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ABSTRACT 


The estimation of parameters from a satellite 
communication system 1s often through the use of Kalman 
memecring. In this work the location of the eye of a 
hurricane is estimated from satellite observations. A 
comparison with a posteriori meteorologist's analysis 
was attempted. An adaptive gating scheme was employed 
in the filter to accommodate "manueuvers" of the storm. 

The observations were at random intervals and also 
came from several different sources (aircraft and land 


based radar as well as satellite). 
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Pee NE ROUDUCTION 


Satellite communications, especially digital, are well 
under way. One of the many observations that satellites 
provide are the meteorological. These give the location of 
a storm or a typhoon in terms of geological coordinates. It 
is possible to have images of the store or typhoon cloud 
cover. 

These observations occur quite randomly in time. 
Meteorologists like to present their forecasts equally 
Spaced in time. 

The attempt, here, is to try to compare an adaptive 
Kalman filtering algorithm to estimate the location of the 
eye of a storm with optimum track values that are given from 
the meteorologist's analysis. The filtering process combines 
all the available measurement data with prior knowledge about 
the system. It produces an estimate of the location of the 
Storm (latitude, longitude) in such a manner that the mean 
Square error iS minimized. The parameters, which are essential 
design elements of a Kalman filter, are the measurement 
noise covariance matrix, R, the excitation covariance, OQ, 
the initialization covariance of error in the filter itself, 
my O), and the transition matrix, ®. 

In most physical processes that one desired nO Meee: 


many of these parameters change during the process of 


tracking. The measurement noise associated with the 
observations can change if a different sensor is used, or 

a Similar sensor obtaining measurements from some different 
geometry relative to the object being tracked. If the object 
being observed is acted upon by external forces, then the QO 
matrix should be changed to account for these external 
excitations. Most processes being tracked will change in their 
dynamic characteristics during the observation time and so the 
transition matrix ideally should also be changed. Further, 
the time between observations quite often occurs randomly 

in time. All of these things bring about the need to change 
the Kalman filter to adapt as the process changes. 

Tt 1s possible to change the parameters of the Kalman 
filter by sensing the error between the observation and the 
prediction from the track of what that observation should be. 
If a gate is established representing 95% of the normal 
random perturbations of the process, then when this error 
exceeds the magnitude of the gate, one can reasonably 
ascertain that the filter is no longer properly representing 
the observed process. In the work attacked here, real data 
was obtained from satellite observations and the qualitative 
observation error was established (PCN #). 

The error covariance matrix in terms of error eliiesoree 
along the track gives a measurement for the worthiness of 


the algorithm parameters. 


Li SADA ere eR iNG TECHNIQUES 


In a linear, discrete system, the state and measurement 


equations are given by 

pe) = Ole) x(k) + Ewilk) 
and 

easy Hik)xtk) 2 vith) 


Pee sewwets the state; 91s the transition matrix; IT is 
SeeeereattOn nOlSe Matrix; H 1s the measurement matrix; w 
icmv~edcee the Cxcication and measurement noise correspondingly, 


assumed uncorrelated, zero mean white Gaussian: 
TC 
Eltw(k)-w (j)] = Os orre 


and 


E{v(k)‘v (3)] R(k)6,- 


J 


and 
Elw(k)] = 0, E{v(k)] = 0 


where Q(k) and R(kK) are covariances of excitation and 
measurement noise. Now if x (k) 1s the estimated state 
value after the kth measurement and x (k|k-1) is the 
predicted value of the state before the kth measurement 


we have: 


The 


and 


The 


and 


The 


and 


The 


A 
x (k|k-1) = (a e where f is any function. 
state error vector is defined to be 


X(k) = ¥(k) - x(k) 


= 


the predicted state error vector is defined to be 


A 
= 


x(k] k-1) = x(k]k-1) - x(k). 


cOvariance of state error matrix is defined to be 
P(k|k) = E(X(k)-X'(k)1, 

the predicted covariance of state error is defined as 
P(k|k-1) = E(&(k|k-1) -&" (k|k-1)1. 

state excitation matrix is defined by 

Q(k) = [(k)Elw(k)+w (k)]-P" (k), 

the measurement noise cOvariance matrix is defined by 
R(k) = Elv(k)-v' (k)]. 

Kalman filter equations are: 

P(k+1[k) = 0(k)P(k|k) 9) (k) + Q(k) 

1 


G(k) = P(k[k-1)H" (k) (H(k) PB (k[k-1)H> (k) +R(k) 17 


P(k|k) = (I-G(k)H(k)]P(k|k-1) 


10 


R(k]k-1) = 0(k)X(k-1|k-1) 


B(k|k-1) = H(k)X (| k-1) 


~ 
-_ 


Six|k) = &(k|k-1) + G(k) [z(k)-2 (k |k-1)] 


|x 


faiemwinittial condition of P (error covariance matrix) and the 
Petesmhe Matrices are the determining factors in the filter 
Pmeeeeeure., For © having main diagonal values greater than 
Emieans that we hayes gqrearer uncertainty in the state 
estimate than in the observation. Thus the new state 
estimate is more dependent upon the new measurement and less 
related to prior estimates. The inverse is also true. For 
Eee@aving greater Giagonal terms indicates that the new 
measurement are subjected to stronger corruptive noises, 

and so should be weighted less by the filter. The gains (G) 
faeemrower, The P is responsible for the initial transient 


performance of the filter. 


eet 


Tif. ERROK ELETESemes 


The error covariance matrix in each stage of a Kalman 
filter process gives insight into the quality Gf tne wteaee 
OCCUrInG. 

The diagonal terms (Pia and Po5) are the variances of 
uncertainty in our knowledge of latitude and longitude. Their 
respective off diagonal terms are covariance between 
latitude and lonertude: 

The square roots of the diagonal terms gives us the 
rms errors in out estimates of longitude and latitude. Having 
the definition of the structure we are dealing with (in our 
case the satellite observations) and its uncertainties 
(expressed by the PCN number-actually by the values "oma 
we can see how the K.F performs through its error covariance 
matrix. Expressing the P matrix in an ellipsoid of constant 
probability, one obtains a visual appreciation for the 
worthiness of the algorithm parameters. The representation 
requires that the errors are normally distributed. 

The joint probability density funcervengus. 


71 /2e" (k| k-1)P > (k|k-1) e (k| k-1) 


where e(k|k-1) is the predicted state error vector. Setting 
the exponent equal to a constant value, we are going 


to have a curve which is an ellipse. This ellipse, 


eZ 


however, does not have its major and minor axis aligned with 
the coordinate system. Instead its axis (x', y'), comes 


from the following transformation: 


Se xcos@ + ysing 


ee ycos@ - xsing 


where x, y are the latitude and longitude in our case with 


Be oe ean Soe 


2 Z 
TL eaeS 
x y 
where cos(x,y) = Pi (k k-1) 
o? = P__,(k{k-1) and o? = Peek | ieas) 
ae de vp eee 


The new variances are calculated by 


Z 2 
2 x as "y cov (xy) 
o' — + ° 
ve Zz Sin2 @ 
2 2 
oO + Oo 
foe x y _ cov (xy) 
ey 2 Sango 


Incorporating the above equations, error ellipses are 


presented in subsequent figures with the satellite tracks. 


i 


IV. SATELLITE TRACKING-SIMULATION RES URe2 


Data from the Annual Tropical Cyclone Data for the 
Typhoon-Nelson appear in Tables 1 and 2 in terms of 
longitude and latitude. The best track data appears in 
six-hour intervals for twelve days and the satellite fixes 
(observed) in random time intervals. 


The data in the K.F algorithm parameters are: 


Ht 0 DT 0 | 
} 
0 1 0 oe 
O— 
0 0 1 0 
0 0 0 0 
pt’ /2 0 
0 DT /2 
[ = 
DT 0 
0 DT 
10° 0 0 0 
0 10° 0 0 
P (wi ; 
0 0 10 0 
0 0 0 107 
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a aOe = 


In Figure 1 we have a representation of the best track 
observed and the K.F track. An error ellipsoid at the 25th 
Stage of the process appears. It seems that K.F follows 
the observed track more closely than the best track. This 
is due to the adaptive gating and Q relative to R. 

The gain history for G(1,1) fluctuates between 0.7-0.8, 
never arriving at a stable value. The G(3,1) reaches a 
Steaole value of 0.12. 

The above appears in Figure 2. The track for prediction 
appears in Figure 3. 

The errors EBl and EOB1] represent (YH-BLAT) and 
(LAT-BLAT) respectively. They, along with EB2 and EOB2 
appear in Table 3 and Figures 4 and 5. 

Implementing Julian Time, we have a comparison in common 
time for 24 points only. It can be seen that in terms 
of the latitude error the K.F is close to the best track 
values. 

Pomsieatecempe, tor better performance the © 
matrix was changed in the algorithm. Also the innovation 


errors, in terms of latitude and longitude, exceeding the 


iS 


magnitude of the gate (,/P(k k-1) + R) resulted in changing 
the values of the gains. The above correction makes the 
filter more adaptive now and the errors EB] and EB2 appear 
smaller on an average. The representation of the above 
correction in terms of trajectories, gains, predz@eti1ons ane 
errors appear in Figures 6, 7, 8, 9, and 10 anGge@iajee 


The computer program appears in Appendix A. 
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LAT 
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3 e30 
Se 80 
5-80 
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BEST TRACK DATA 
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2200-00 
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2212-00 
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2500.00 
2506-4 00 
25122e909 
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2606-400 
2612400 
2618-00 
2700200 
2706200 
2712.00 
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2838135200 
2900-00 
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118250 
117.80 
117-20 





Foe ed 


ES 


XH) ; 
ack values (BLAT, 


pier. ee, 


(YH 
PSt tr 





Koa 
B 







ye aa 


128-1 13258 “1a725. 1424s 8 Seems enol 
LONG — XH 


Zee 


Trajectory of Storm Nelson 
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ERROR VALUES IN TERMS OF LATITUDE AND Leer: 


ATT “a= =A FRO ey FNAD 
FINKRGAD —-7.799 0 «nA —) ef) ats) 
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Trajectory of Storm Nelson 


Figure 6 
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V. RANDOM TRACKING 


To examine the adaptability of the K.F, another storm was 


created. The equations for this simulated storm were: 


BLAT(s) = BLAT(s-1)-v.°T 


BLONG(s) = EOC ERD) ad 


LAT(s) = BLAT(s) + V(s) 


LONG (s) = BLONG(s) + V(s) 


where T=6hr, Vi =10°/24hr, eo and V = measurement 
noise (created by a random generator subroutine). 

Implementing the above equations, a "true" and an 
"observed" track were created. These two data files were 
named BESTRACK and OBSERVED and appear in Table 5 and 6. 

Simulating with the above new data, the K.F algorithm 
appeared to track well. Figure 11 shows the "true", 
rObserved" and K.F tracks. 

Figure 12 indicates the gain history in terms of Gll 
and G13. This approach gave stable values of 0.27 and 
0.01 respectively after the 4th discrete point in time. 
After this time the gain does not vary any more. This 


means that the innovation error, (2 (k) -& (k| k-1)) is 


3} Il 


weighted each time by the same quantity after the 4th 
observation. Having lower values in the diagronal terms of 
the Q matrix in comparison with R, in this case, means that 
we have greater uncertainty in the measurements (observed 
data) relative to the model uncertainties. So the gains are 
smaller and the filter no longer "tracks" the measurements 
closely. 

As far as the latitude and longitude errors are concerned 
it can be seen that EB1l (YH-BLAT) and EB2(XH-BLONG) are 
nearly zero and are smaller in comparison with EOB1 (LAT-BLAT), 
EOB2 (LONG-BLONG). This shows the ability of the algorithm 
to follow the "true" values more than the measurement if the 
latter has been corrupted with noise. 

The error values and the plots appear in Table 7 and 


Figures 13 and 14. 
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BESTRACK DATA 
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OBSERVED DATA 


LAT 
160200 
157.230 
154.92 
1952.71 
150.023 
147e83 
14495 
142.51 
{39.99 
137.26 
134.62 
132239 
130e31 
L277 aor, 
125-19 
122.44 
120.219 
117e40 
114-87 
112.92 
110.03 
107.65 
104 e83 
102.96 

99.90 

97285 

95 e29 
92.44 

90.51 

87218 

B84 e834 

= e40 

30.00 

77278 

75 208 

72250 

69 «96 

G67 waz 

64.79 
G27 5 
59.55 

5S 7 aa 7 

52.43 

50.06 
47291 

44.38 

42.53 

40.14 

37 e2b 


34 


LONG 


32830 
2260 
1-72 
21201 
—-OQOel7 
—-1.Q7 
—2e25 
—-3.09 
—~4e.21 
—5, 44 
~6.58 
—7e31 
~—~7.239 
—-9.33 
—-10.01 
=l12a0 7 
—-13.30 
—~14.233 
—“1606e17 
-17.05 
—~18237 
—-183.7% 
—-20.-30 
—2leoS 
—-21.91 
—~23-258 
-23.59 
= 22 
—~262e36 
—-27.30 
—-28.20 
~28-92 
~30el1l2 
—~31.20 
~ 32.24 
-~33.%8 
—-34.41 
—~34.95 
—~ 36037 
—~37.213 
—- 33.743 
—-39 ..27 
—-40e1+4 
-~40.79 
—~42e32 
—-43e212 
—-4%.06 
—~45 244 


ue ry ri a ee re rt eee ett ae > te ey tee wee ee er ae er ee ere gee Se ES es ~ 


, LONG) 






XH) 
t track values 
BILONG) 






me K. 
e Be 


@ 


Fars 
_ 
vel 
in 
. 
ed 
=F. 
rena 
rake 
ne 
B 
Ld 
E 
4 
Aap 
Bi 


091 £°Zo1 b°SE1 1°E21 BOLI S*86 2°98 BSZ 919 L*6rO"ZE 
| Bec S| 





55 


ae 


a0 


= 2 aE eee We cal] 10) ee) 10) Oa 
LONG - XH 


more) 


Solel 


Soy. = aul a) 


Trajectory Of Eretitious. Storm 


Figure 11 


— ee ee lt 


ARTOOTSA pue UOTFTSOgd JO SuTeD ZT sanbty 


AWIT 
iS ea c OF me Sh, Wea = aa iG a ae /, al 8° OT 6S Gat 
WOE NONI EN EN 7 ONG 3 WENN EN Wi, Neal KX KXDSOKC KK Se YZ Eons See aN a Sé nya 2 
» Moo 
pe / 
\S j 


Eee Aa TOO TSA (7 (6) 5: x=x- 
uTeb uoTRATSOd (T‘T)9D?¥x 


870 £0 9°0 S°0 ~7 ORs Gite 
[oe 


60 


af 


36 
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Vibe CoNChUSTONS 


The K.F approach to estimate the storm's location appears 
to be very accurate. That comes from the comparison with 
the meteorologist's analysis results. Concerning also the 
fact that the latter was performed after the storm's occurence 
one can see the advantages of the K.F algorithm. During 
operation of the filter the actual residual sequence 
[2 (k)-2(k| k-1) ] is compared to the gate [VP (k| k-1) +R(k)] 
which actually is the square root of the residual covariance. 
Being a white Gaussian sequence, the residual is bounded by 
this gate. 

When the gate is exceeded the model is determined invalid 
within the filter and a modification in terms of Q and G 
takes place to adapt to the situation. At this point, and if 
the excess occurs in only one component of the vector residual 
process, one can further deduce that the measuring device 
generating the particular component is the source of 
difficulty (a sensor failure). 

The error ellipsoids of the process also give insight 
into the filter performance in a more general case, referring 
to many sensors with a greater variety of uncertainties, the 


adaptive K.F algorithm could be a very advantageous approach. 
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ie sins) pe homeda)yoCn (J) ,J=1, 42) 
2) 


niet ba Wk, bo. 0) 


Peeve See Oy Pe NEXT IRR SE LINES TO ELIMINATE 


Pee ana CEhkCK 


Ae eA eK ok ok 


meee en ela 1), Gad (1)),LON> (1), ELIAE (T) 


Ter (0), 
: bree Nickie.) .2C 


(I), I=1,NZ) 


eee ye leeenn i pT (1), LONG (1), PCN (I) ,I=1,NZ) 
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CY ie ae 


©) 


Clie UGA? BE “00 


Cy 


) 


~ 


O 


Co; Te i, ees 


Corel! Cy 


Ge 26- (Ge 


PGi. > 


- 


1 
* 


s/i4 


Yl, 


IT=k 


BLCONG (1), SIND (1), [=a 
WRITE (9,374) (GZ IME(1) ,LAT(1) , LONa(ie eee 
Pena (eee 
BEN (1) 70 =a 
FORMAT (2% perme) 


THiS PROGRAM COMPUTES TRE FOLLOSING RAD 


G(K) = 2(K/*=1) 4at* (Hee (Kea) ee 


ae 7x) (I= (Ky) area 


J) DEFINES THE COVAR TANCE OF TAs PEA Sane eee 
EACITTAITLOCW OF “Die hee 5 


J) DEFINES THE RANDOM (GAUSSIAN) MEASJRZAENT NOISE 
{HICH IS ADDED IC INE C8SERVACLE SIGNALS 
q_(1,d) IS THE FORM pay ge ae 
TAE DISCRETE POINT IN TLiME,THE STAGE UF THE PRCCme 
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C5 


c Pewee eee ei C OV FEE On AT A wIVEN A SAZELES) 
S 

e 

ee ee ee ea, (COV ERROR AT X GIi¥V2! K-1 
oat et eS) 


e- eterm oon, go = numee* CF CCL. ,CEBS. 


rz 
rrt 
f+ 
od 
| 
(74 
re) 


c 
. 
c Wenders n Gee ti brRAZTLUNS Cr 1H 
S 
C 


C FR RRKKEKK NEY NE YS 3 Xt 58 OK a a te KR RK RR OR OK RO KO KK KK KR OK OK 
WET = (8,717 1) 
yi? ft PORSAT CO/ ok TK OME AR ACK OR KARE KK TY EW Nae a ee ee 
#OX, 1THE STORA PRACK INPUI VALUES ARE AVLALABLE',, 
HES2 AkE ECHC ERINTED IC THE TERAINAL AZ TEES, 
ee TSENG 7, 


eee eer Neem ENDS NCPMALLY AN InvUl FILE 


1 


~ 


; 


© 
U 


Age 
Cow stoi Ue LHe raOGRAY =—= ALS 
X 


i] 


eee 

#2X,'BE PRCLUCED THAL AAY BE PRINIECD G2 USZD 

1 FGR INZUI',, 

2X, "SING THE SAME FORMAT STATEYENTS TO READ as 

1 WERE), 

*2X,'USED TC WEITE CN UNIT 4. S23 WITH a EEX 

Me oe DenI TONS "7, 

#24,'THE BRANCH AFCUND THE INDUT CAN B2 US#D. 

1 PRESENILY',, 

42X,'A EILE CF THE INPUI DATA IS BEING PRODUCED AND 
(cert, 

*2X,'BE FOUND AS --K CUIBUT A-- CN YOUR A - [uSKk.") 


cay 


C 
Pine ieeeoin s OUEST EON TC THE TESSINAL --ASKS IF 
S iMG wier tbe eiS 10 EE USED 

C 


44 


GC Cx. 663 


oO 


1235 FORMAT (' KHEBE CTHERWISE S 
1916 FORMAT (140,5X,5H X = ,13,/6%, SHGAINS) 
1992 SREDE(d, 10) 


2X, 'OG NCT USE Al MEY TC EAL eae on eee 
1 YCU ea ee 

¥24'*TO LOSE YCUR INPUL Fite THAD 2s oe 
LUNs, 

Aw RK AK KR RR KK AK KAKKR SAKE KAS ES KK 
*2X'FOR CHECK CUT ACTIVATE THE GO IG B@ogbeSIATea ae 
1) “Cx 

*2K' CINE 135 LO°2ND {teh ult eek ee 

MOU RR RR RRR EE RRR ER RR RRA RRR RR RRR RRR 


*2K,'DO YOU wANT TO DEVELOP THE INPUL Pate eee 


(3,33) DAWs 
i 


33 FCRMAT (A4) 


7717 FCW (181) 


W El Seo 2) 


1236 FORYAT ('ALL INPUT SHCULD Be FLICATLAG POR ee 


WRILE (3, 1235) 
DECTE UE stews 


10 FO24AT(/,5X%,' THE DISCRETE KAL4AN OE LO eee 
WRITE (8, 30) 

30 FORAAT (5X,! ENLES THE OSDER OF THE SYSTEY 
ee Cmca 
4EAD(S,40) KN 

uQ FORMAT (11) 


45 


bale 

oe 

orto L= 15% 
ary =o 

Fema dced- 1,1 
Heme) = 0 


poeta 1 
Bee (ied) =O 
Em (ioe lye= 4 < 
Poli) =O 
Ene = 
fae A (I,J) =0. 
WRITE (8,410) 
WRITE (8,7771) 
7771 FORMAT (5X,'CO YOU WANT IC COMPUTE PHI 
7 LINE FRCM!, 
Oo es 
READ (5,7772) TAN 
7772 FORMAT (A3) 
Me PAs EO. 12H) GO TC 7773 
TF{IAN .E.TEN)GO Tc 15 
Moga .RLTE(8,7774) 
ga POuMAT(//,5%,'=NTER THE A MATRIX") 
DO 7775 I=1,%8 
DC 7775 J=1,% 
MPITU (a, 7776)1,0 
io pee {54,'A(",21,",',11,")=") 
Pan Deore 76) Al) 
7778 EORMAT(F10.3) 
7775 CONTINUS 
WRITE(9,7779) 
WRITE (98,7779) 
feo FOeqAT(52,° THE A AATRIX ) 
wor eeu T= 1.N 


«4 
G3 
wT 
texe 
Ls 
a» 
Cc 


46 


1734 


T78S 


70 


Aas (57 SON eae 
PO) ee ye, 
44772) 


j=1 
Wel Te feo yeh. 
WELTE (4, 
TOSWaT (55) ee vas 


Tit Clits io. 


mr 
(44 
oa 
eho GD 
Ww 
~~ 
oS 
— ft 
Nh 
dal 
a” 
rea 


riy 
©) 
ry 
ae 

~} 
b~4 


eee” 


tr 
— 


ical 
Pee AT (5x, ' Ge 
DO 77382 1=1,0 
POU) o2 c= 
PRS (oe vos 
FORMATS <q! Ga ee eae 
READ (5S, 77 Peale 
CONTINUE 
wBITZ(S,77E4) 
WSLTT2 (8,774) 
POMVIA DT (5c; 

pc 7785 
WRITE(9,90) (5 
WHITE(3,30) (E(I,J), 
COm Caco 

WRITE (8,30) 
FGgelal (6%, "ENTER 


ply ane 


Lee 


i = 
be 
te 
mele 
4 
os 


m) 
ela 

Te ee 

l= eee 


LA Se ee Nc meee LkhaAtviol ou 


MADRIX--PEI') 
PO ie= 1 
pO 1 J=1,N 


WRITE(8,00)1,J 

FOWMAT (5X,'EHI(',11,°, *ol eee 
BE (ou oy Orr Hiei en 

FUR AAT (F 10.9) 


All 


| Cohen: 
foe te, 30) 
CeCe) Ck, Lee eA LT RIS (TRAN 
Deg [= 1,8 
Pa etme mere One eh te ed), J=1,N) 
So FURMNe2T{1P/e11.4) 
eee a le 7 el 1. 4) 


Coe ok ee oR RO Rd OR OK OK K 


tw 
W 
Hi 
ti 
M4 
Cc 
i 
ar 
4 
je) 
in| 
_ 
fe” 


GO TO &88 
CK AR ROK A KKK OOK Ke OR OK kk KOK OE OK A 
BOUO SF FITE (eG, 100) 
felt Cn Ato een CUNaaw.. LO CAANGE ANVOPELENENT Cr LAE 
lei Tht) 
our eA 45,1719) LAN 
110 FCEMAT(A3) 
Pies eee TZ)GCTO. 13 
Pea etteerr t) GOLTC 1000 
WeLi LE (3S, ao 
fee OE AAT (DX Peso OF SS SALE LS. JU yOu SANT 
1 TC neg: 
Peers few RE FT 1S Lae kCOW AND J IS Iz 
aCe wor N ce) 
Br (ogo) Lo 
foo FORMAT (2m 1) 
Nee ee (d, 50) 1,J 
READ (5,70) PHI (I,J) 
eile (oo) 
DO 3 L=1,N 
Meee) olen b{ ly) ,J= 1, ) 
WRITE (d, 140) 
fer wPOITAT (>X,°ANY CTHER CHANGES?') 
Geter 3 U0 
howe (97.141) 
ooo rere te ek, 'DISCAETE TINS, 


Gre) 
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Wit 2 pono ee 


1117 “ARETE (9, 91) (See ey 
WRITE (8,410) 
WRITE (8,150) 


150 FORMAT (5X, '2zNTER LHe DIAEWSITON CE BAM Deri oeee 


rej 


1 VeSECR (Nee 
EEAD(5,40)1 
Mn tO. oO} 
TOO FORA (500 ee 
+ Ke ee kee eee 
po 4 I=1,N 
BO Wood =A ek 
PETTERS lO iss 
170 FORMAD (5X, 'GAYMA(',I1,',',11,') =!) 
Nae TO aaron 
U- CONTINUE 
Ri eee eC) 
toO FCRSAT ('0',5k, 'THE GANAA MATSEY (Distalaeee 
1 MALRIX)') 
DOs ow l= el 
5 wRITE(8,99) (CEL(I,J) ,J=1,L) 
V1 Oe wee TE (oy 0G) 
910 READ (5,119) TA 
if (IAN.2O.12Z)GoTC 29 
IF (IAN.NE.1Y)GOTO 1109 
melts (354126) 
Beenie (O70 oes 
errs (8/170) 1.0 
RE AWS, 70) LER Cleo) 


OF THE JisSiklc Jao 


ta 
to 
rH 
LS 
tx 
le 
a) 
& 
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MaATE (8, 180) 
emo ha, h 
paeeeeCce On (bel (143) ,J=1,1L) 
hom Be, (40) 
GOTO 910 
29 WRITE(9, 180) 
Doe igen =a 0N 
imei "Weeer(o, 91) (DEL (1,7) ;J=1,1) 
Coen 1s, 10) 
9910 “RITE(3, 399) 
390 FORAAT (S54, ‘ENTER THE ELEMENTS OF THF CCV OF %') 
C 
MP SG GeusthATueD 19 (1,1) 2 <K 1719 6 4459) 
ee =u Oe 4 


fem eo -2 1a) 
DO 21 J=1,L 
Woelis (Gee 101 pd 
comentario, sett si dt, ',*,i1,')=*) 
READS, 7 ORMW I, J) 
Fee Gi TN UE 


W 
Secret) | ok, tHE COY OF 8!) 


Pee WRiDS, 20) (8 (1,3) ,J=1, L) 

1400 ®RITE(8, 100) 

940 RPEAD(5,110)IAN 
PeeNNe Boel Z) GOTO 89 
IF (IAN.NE.TY)GOTC 1400 
KEITE(8, 120) 
READ (5,130)1,d 
MERE Ce, 310) 
Reo TONG (1, J) 
ARITE(¢,321) 


0 


Peo, 7 ee CF ie) 


ew, 
Ca 
> 


1115 wRITZ(9,91) (W{I,J) ,J=1,L) 


1395 FORBAT (5%, ENTER TBE ELEAENIS Of DHE (Gee 


PCn? aed 
Woe (di 0 0 )led 
290 SCR4AT 1S, * Hee ete tee 
PEAD(S,70) H(I,d) 
TCO Nie nz 


Aye (go ee 
210 FORMAT ('O',5X,'THE H 4ATSLX (OSSERVATICH MATAIX)*) 
BC 2 We hep 
Sete (oo 0} (Ht lide ee 
ez COON Ree (ee 
Or eee to eh 
PP(lane ees tence ta sc 
IF(IAN.NE.1Y)GOTO 1200 
Oe 2 (ye GC) 
READ (5,130)1,d 
WRITS (8,290) 1,0 
PE ao 9 0.) aati!) 
WRELEXSG | 240) 
BO 3.) 1-48 
G Hele (gn90j) 7 (ae eet 
WRITE (8, 140) 


mal 


Giga 9 20 
SOM kee (Ss, 219) 
Gta. C= et 
eee st oo, Ole () = 1, N) 
WRITE (3,143) 
Wee (eee 70) 
270 FORMAT (5X,'ENIEP LHE ELEMENTS OF THE MEASUREYENT 
WewO Tseity 
#14,"COVARLANCE MATHIX--B ') 
eee lo. f- 1,4 
Be mel oma 1a 
feats, 2 e011, wl 
Pee neta, tn (', 1 1,', pe t,’) =") 
READ (5,70) &(1,J) 
116 CONTINUE 
pete, 200) 
Mo LEE (9,290) 
POO EORTAT (4,52, '2EHE MATRIX {(AEASUREMENT', 
44X,"NCIS= CCVAPLAHCE MATRIX) ") 
femecs 17 1=1, 6 
ee (eee (1) oe J= 174) 
fev write (6, 109) 
gan) READ (5,110) 1AN 
Riise GOLG 79 
Pinney GOnC. 1.400 
ARITE (8,120) 
Beemer o) TJ 
Pieter 25 0).2 50) 
ames O) nd lJ) 
WRITS (3,29C) 
SC wise f= 1! 
oma atne (8,90) (R(i,J0),J=1,% 
wRITE (8,140) 
GOTO 930 
Oe. WA selena 


rr) 


ow 


1114 


da 
Cas 
ae? 


se 
N\.? 


51 
44 
Je 

150V 

0 


36 


21 


PELs) 


Zuo 


+7 Glo) 


ao ee) Cate 
WRITE (3,919) 


gi) p= lie 


WELT=(€,2e0) 

FORAAT (/7, 52, Utes ERX ant ee ') 
PCueen i — (ee 

pe Sh J=17 8 

ARITE(8,34C)1I,J 


FORMAT (SX, "EXXH1(', 
READS) ONS a eee 

CONTINUE 

WRITE (8, 351) 

WRITE (9,351) 


pg. eee 


PCR4AT (0, 5%) Pere 17 i) 
BRITE (0,99) (oReee yd) p= ee 


JRLPE (8, 100) 


Reap (5,110) am 
LECTANS Ge 2 2pee CO. 
IF(IAN.NE.L¥)GOTO 1509 


WEEE Ge, 120) 

PEMD (og a) ad 

Se Ere yes) 

BA Deo ge) oe el Cle} 

reas joel) 

pO 359. f= 15.0 

“Ee ULE (3,30) (2844111 ee 
we Lats 7 14) 

GORO 3950 

peo 11 Vo) teal 

WRITE (9,91) (ES AMR o od) = ley 
Re CIR (95 43} 

nos (Ola o ) 

FORMAT (94, "ENTER OTH Ee ae 
Phe GR iD as 


r4 
iy 
LO 
ret 


ete Ae eee 


VALUE'} 


5S 


1,313 


V665 


511 


m™ +) 
5 


Ui 


2019 


Pee o 7) 
Perth. (1 2} 


iT = (fide (1)-1400) * (201) 
IS= (LINE (Z)-1400) * (.91) 
T=2UeIT 


S=24*IS 
T=TIME (1) -IT*100-140041 
S=TIAE (2) -I1S#1)0-1400+S 
Mt=sot 


CONTINUE 
Jh=1 
HPITE (38,7665) 0% 


aU Creer, (4 } 


WRITE (8,511)N,M,L,ND,45,LD,UN,DT 
FORMaT (2X, 2HN=,15,54,204=,15,5%,2HL=, 
ees Gon. T.32 
#5X,30Lf=, 
poe > CG ckwt is 1 Soon GHET=,6 10.4) 
WRITE (8,532) 
FORMAT (/! MATRIX R ') 


Dome el =e 

emt e(s, oO) (a (led) ps=1,%} 
WRITE (8,544) 

FORMAT (/? MATRIX @ ') 
Heese = 104 

WRITE (3,90) (0(2,J),J=1,%9 
ARITE (3,555 

FORMAT (/! MATRIX PKK41") 
DO 30619 I=1,% 

WRITE (3,90) (QXKM1(L,J) ,J=1,N) 
[F(IANS.NE.IPH)GO TO 6789 
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Poy ooe iN DS 


O77 59 


r 


€00 


3920 


44 


21 


S25 


444 


(Tia ( i ee ee oy 
Mee (Lida) 1 eee 
ol 


Pe 

S= 7s 
TE=TIME (1) - 114 105-1 
S=Tiae (2) —-1S4 100-4072: 
pr=S=71 

JA=2 

Welle (G47 a0) A 
CONE u ey 


Get h (5, 7878) 0 
FOR4AT (5X,£10. 3) 
CALL PHIDEL(DT,N,L,A,2,Pul, DEb,01,02) 9 een 
WRITE (d,€66) 
FOFMAT(/! a ') 
DO eo) la 
WRITE (O60) (Eur a 


PORMAI(/! DEL ') 


we ite(s, 90). Dib D aty een 
CALL TRANS (DEL, A, 1, Dent, De 
CALL PRCD(DEL,DELI,N, 1,950 ) coe 
CALL CUNST (WUT, 1) 0, de oon 
HRITE (8,543) 
DC J025) TaN 
WRITE (3,90) 10 (1,0 l= en) 
WRITE (3,054) 
FORMAT (/! i ) 
DONS Ozen I= 153 
PPTPE (cd, 90) (H(i 3) ee 
Jie 3 
ae mie verona i 


a 


co 


Ba27 


Bead 


a2 3 


eae 1 ip ll, 1D, LD, SveDe 
Boe ar? 240, 2 1G. 4) 

BO 5327 (i= 1,.N 

peu) tA (l,J) ,J=1,8) 

pom oo7s I=aeN 

meer e (9,09) 6 (2, 1) 

ve 5829 ease 

met Pee Donate ao), J= 1, N) 


Doyo 30 b=, Ss 
Poe Oe) ean 1 (2,0) ,0> 1, N) 
Demon el) i= 
eee ee One (HC 1,5) ,J=1, il) 
VORS 72] == i, N 
Aoi of 
DeRose = 1, 
See ow). (eetl yd) ,d=1 70 
WRITE (3,7777) 
Oi ee oe 


3 (a OP eee ee 


eee CSC ee ok OTN OCR XXX ' 04) 


Wee cae 7 eG ord K 

Peat od dl) N,%,.,ND,4MD,LO ON, DI 
JK=HIG 

eee een) 65.5 nih 

Pie ol ieee ©, NE ,ME,LD,NN, CE 
JKx=4 

iRITE(8,73c%) JK 

Poe Poked = 17 is 

peat die (AT ied) pJ= 1.8) 

JK=432 1 


-—t 


© 


WRITE (38,7865) JK 
LO 7236 [=1,N 
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7250 BEAD (4,91) Ea 
DG 72S seal 

Tz23d FEAD(3,91) (0 eee 
DO 7233) Pain 


READ(4,91) (PREM U (2 goede 
7239 @RITZ (8,90) (PXKAT tt pee c= 1 
De 32h T= 
7326 BEAD(H,S1) (iia eee 
PO a2 SE = 
7329 READ (4,94). ae 
Ow 32 0) = ies 
7327 Ree 91) (2 eee 
iT= (TIME (1) -1300) * (.01) 
is= (1aa8" (2) ~150G) = 04) 
T=24*1T 


=TIAE (1) -IT#100-1590+1 
S=TIMB (2) -IS#109-1500+S 
‘ ; 


rip eso 
JK=5 
RITE(38, 7065) JK 

CALL TRANS (L EGRSEIy 1 eet ey) 

CALL PROD (DEL ALELT 1, Nee 

CALL CORST (841, 1) 407377 2D 

CALL PHLIDEL(DI,%,L,A,B,EH1L, DElgee 1, ol, Nnen eee 
7234 CONTINUE 

p= or 

CRITE (3,91E&2) 
9182 FURMAT (5X, ' ¥CU ABZ AT DO 2222 ' } 

C DLDLDLDLLDDLDDIOLELDLOLIOISI ULL Lome 
BB=1500+IT#190 


KKK=1 


ae 


o 
Cc 
eS 
= 
1930 
133 3 
1984 
i265 
ie 
C 
C 


, ARB: 


Cy 


CY AIG ©) 


Pee CATES SUEEOUTENE TO CALCULMTE JOLIN PME 

eke FOP EVERY SiLG84% POSIII 

LLLMESINI (TiM? (K)) 

IDAY=ITI41£/100 

THOUn=ITIME- (IDAY*109) 

LF (IHOUR.2C. 90) IHCUR=24 

MCD6=S40D (LHCUR, 5) 

CALL JOTIME (ITIME,JULEF) 

LF (MOD6.EQ.0)GOTO 1986 

MEIMB=DNT (EE) 

SDAL=MIIME/100 

ADAY=4DAYV#100 

MHOURSA4TIME-MDAY 

IF (MHOUR. 32.24) MPIME=ITIME 

CALL JUTIME (MTIME,MJULHR) 

CO sueenielo 

ARITE (8, 1989) JULUR,TTIME 

HOM 7 7° GUnTNN FOUL 1S *,19,' , ACTUAL Tisé 
>) 

GRITE (8, 1984) 4TI4E, MJULAR 

FORMAT (' MCDULA 6 TIME= ',15,"' ,COBZESPUNDING 


femurs p= 6 1S) 


CG ans 


Vero LAN LEME nCUlLlnNe 


<7 GME T ADE SODULA 6 JULLAN TIAE FOR EFKE,7KM,XKE C 


*** AQDTULA o FOP JNLIAN TIME 
LF (MHCUP.EC.18)GOTC 3187 


IKM72(K) =4JULHE-72 
TKM48 (KX) =4JULHR-48 


a 


Cx ViGiP 263 


3135 IF (PCN (X)-2.NE.9) 


TXM24(K) =1IDlea—24 
TKM 14 (K) =4IULA§-12 
TK 46 (K) =40CLHR-6 

[X26 (K) ig dbakt+e 

TKE12 (K) =40 Utes 
Te 224g (ky = JUL eee 
IKPG9(K) =“MJULHE+48 
TK272(K) =AJULER+72 


167 SCG ide: 


45% OND JULIAN ©] Ee ee ee 


i? (PCN {KE )— 5. Ocoee So 
aly cre 
Rote =e cs 
GOm ors 154 


313537 LE (PCE (Kk) =—SaNEl0) (Gem ee 


Ba pelea) ces) 
Rey c= ee ae 
Ge O34 


Ga 
C) 
ry 
Cy 
en 

ened 
la 
£ 


R(t, 1) 2.03 12 
Rae = 2s li 


S130 COs 


jee (BS 5) = eee) 


5445 FORTAT (/////7 20k, eee 


Ho, S895) re sees 


9313 FORMAT (3%, 'R=',13,5X%,"HER=', £4. 2,5%) eee 


dee dee re) 
RPITE(4,9313)%, 50,501 a eee 
WRITE (8,313) PCN (2), Daas len 
313 FORMAT(3X,' ®CN(X)=',F10.4, 34%, "£1 = "eee 
pre) 
HRITE(4,312) PCN (4%), D1, (1,1) 
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Be 


oO 
~ 
LO 


ey 


G50 


0 23 


3819 


as) = oe 
erie oss See att, 1} 
a ee e/a 2) 


see) = oe ee at = 2} 
Sey = Dix=24u 02, 2) 
WRITE (9,795) 
FORiaT (75 MATRIX C¢ ') 
Dew oes 
ura meld) 7 J= 1, N} 
pale ye Cero cr 
fe pee 900 
Gene md’, i, ,A,8,2H1,DEL,01,02,N9,49,L0)) 
JRITE (8,975) 
FORMAT (/' PH ') 
WRITE (4,578) 
Se ngo 7 ml 
Pome ee (sel (i,J) ,J=1,8) 
BRITE (8,99) (ERI (L,J) ,J=1,98) 
Dey Umea , NN A1, 3,0, 2ul,a,a,",G/ni,00,A40 
,L0D,8) 
WRITE (4,658) 
WRITE (3, €56) 
FORMAT (/! PKK ') 
Demoees i= 1 ,N 
ope, $0) (PKK (1,3) ,J3=1,4) 
0) (PK (ipl) ,J=1, 2) 
Orit ape eet alana al 
CAL BROCE (H,AKK41,N,%,1,2K41,N0, 80,4) 
@RITE (3, €810) 
FORMAT (/! Bana ty 
HPEES (3,90) (ZAX@id}) ,g=1,2) 
eect ueea NA 1 (J) ,J=1, 4) 
TRITE (3,8619} 
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3319 PCSAAT(/ EA eee cee 
qd we lTE (3,50) DAL (eee era) 
At 1) =) ee) 
o( 2) = Oss) 
7 oe Eee OMe) 
WRIT? (4,90) (2) 7 d= 1s 
WRITE (3,8811) 
WRITE (3, 90 (ae alee 
CALL SUB (2, 2K 3 yee ee 
WRITE (38,5445) TINE (K) 

C jel flo pc ome” 
WALLS (83,3020 
POLtheiy, * 
WeTT ECG, SO) os (ld = 1 

LE (he beer) ~COm comes a 
Gn R= (2h wat 1 (2 2 eee lee ees 
Tf (ABS (E (2))-GETE -LI:0.) GOUT eee 

G(2,2)=0. 5% (1524642 a 

d(2, 2) =WoOCOSSh (232) 

oS 229.54 (9.33546 (eee) 

EKK41(1,1) =Z*EKKA1 (1, 1) 

PKAM1(2,2) =2*EKKM1 (2,2) 

PKKM1(3,3) =Z*2EKKM1 (3,3) 

pi (454) =2* Paes) 
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KOK RE OR KK OR OK KE ROK KE OK ") 


Cie, ye ey 


Cen 


Wp TP Vc 2 Gok een en 
W2 ETE (eG, 9 192) GA et 2) 

2203 CCRZINUS 
GATE= (PA Nd 1) lee te eee 

IF (ABS (2 (1) ) -GATE <L2.02) "GCy lem 

Gil, 1)=0. 541.246 (12 
Ol, 1) SIO 00 ss 7 ie 
G{371)=055* (02 3s2cerene 
2KKA1(1, 1) =2*EKK41(1, 1) 
PKS 4W(2,2) —2 Fo Ree 
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Cy (eae ee) 


Ci 


9192 FORMAT (9X, 


X41 (3,3) =2*PKa41(3,3) 
KAT(4 SS 
mCi oe Mec eyea |) 
ie soto) GAL, (1) 
4 


fe gs, ee CemGmeGATE. GagtE= %,F10.4,9%, "= (1) 


pes Ce kK) 


ts 
FIOL4, XXX") 


Cenieeoi= co tn 
Sai C3 51) 
De 22-3 = 1K 


WRIT! (8,90) (2XKH1(1,0) ,J=1,N) 


2A 

cK 
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9191 EOP 
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pe 
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OY 
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oo i 1 


ITE(4,90) (EXKA1(I,J) ,J=1,58) 
Hoe ee oS ) 

WRITE (3, $9) 

ee ey © CAmerels 4 ~) 

OS Be) Ae 

WRITE (8,9C) (G{1,J),J=1,4 
Te) 
Paeeietoyey) ( Zan b{J),J=1, 
HRITE(3,8811) 


fee (4 ae 1 1) 

FORMAT (/! Z ") 

ey ?°) Gnd, clo tt) 
ees ee 0) (2) eda 12) 


ete sOn (Gath, 1,62 ,ND,4D, LD) 
fees S00) (GE), J=1,N) 
Crema DD XK hel, Ge,s, 1,kKK,ND, 19) 


Sameer (Pa. Yak NN, 1,%KKM1,50,40, 


ARITE(&, 6011) 
WPITE(4,6011) 
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PaGemSGADS. GATE= §)F 10. 


LD) 


a4 


3o25 


Jo85 


SS ie) 


oO 


FURMAT (/! XK ") 
PRITE (3, 50) (Xo ee 
wEITE(G, 90) (26h Oo 
eae = aia ae 
Gens) = Ss eee 
Pd iia) =a 
2K 19 (Kk) = 26 cole 
VP (KO) = ce tale 
PEN) = ok ne 
KTIME=INT (TIME (K)) 
RX TIME=FLCAT (KTIYE) 
LF (RATIME.WE.TIAE (Z)) GOTO 8813 
KJULHK (KX) =JULSE 
BOR 1 (hj e1 eh) — Paw) 
EOB2 (K) =XH (K) ~LONG (%) 
CONTINUE 
WRITE (4, 2812) 
WRITE (8, 2812) 
FORMAT (/ XKKM1 ') 
WRITE(8,90) (X484915),J=1,5) 
ERI TE(S, FO (eal eee tee 
THOR = eK Ty) 
KHP (K) =XKK41(Z) 
BB=1500+ 114100 
pO $684 t=1,4 
LE (TIMF(K)-EB.LE.9) GC IC 9937 
EE=3R+6 
IF (BB-TIME(K).LT.0) GC IC 3660 
IF (BE-TI4E ({X).GT.0)GO0 TC 9699 
YKMo (X¥K) =4KK (1) 


KAM6 (KKK) =XKK (2) 
GO TO 2698 
CC=BbB- (1590411 #100) 


Weegee (41,5050) CC, ed 
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pe ooeemeree Ji (oe CC =' rac 4, ' Poses (Uist, * 
Ve eee OD 


=e re 
ies Ys) 
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Pleo oe 8G 3094 
(PUMe tae ieee. 2t. 3) GC TO 9634 
9700 BB=8B-24+1C0 
9694 IF (TIME (%+1)-BB.LE.0) GC IO Sou4 
If (3B-II“E(K).£C.0)GC IC 9698 
DDL=BB-T IME (x) 


SPITE) eomleleca, DDI, IT. 2 
Cerro (aes) bo, 0Dleme, x 
Porc ey, Fe= 9",F10.4,' DOt= ",FiD.4,! Po pee 
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WRITE(8, 4312) 
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WRITE (49,4812) 
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YKMG (KKK) =XKKAG (1) 

KKAG (KKK) =XKK46 (2) 


Jey XKKK=KKK 
Th = 
WRITE (4,8496€S) 

3969 FORAAT(/,' YX46 KKM KK RB 
| ieee 


WRITE (4/90) YE46 (XKK) ,XKM65(EKK) ,XKKE, 85,71 


peat 
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mo (So ei | 
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n (dao 12) 
HRITE (4,90) (44886 (3) ,J=1,8)) cee 
3222 Ll= (2. 4E4 441) = 9500 eee 

LS=({1Ti 42 (at2) — 1500) ieee 


T=25*1IT 
c 

5S= 2075 

L=TPys (XK+1)-1T © perros ee 

S= Pies (A+2)-15* 0G] 1200 se 

Di=o24 

JA=6 
u Q0Ttse (3,766€5)0% 

ARETE (4, 55-2) 
Were eco ers 

c DO 3922 f= 
e ALE (3,50) Cee Ty 7 J = ie) 
Cou MR ITE (4, 90) Geeeeee to =e) 

CALL PRIDEL (DT,N,L,4,E,EH1, DEL, 1702, eee eee 
C CALL. [Re vS URE eG etl cee uN ge oa 
C CALL PROD(CEL,DELI,N,1,5,20 4 Noe ee 
g CALL “CONS 1 Ami, 1) 4.0 3GN; Nt hee 
C LOLLOLDLELOLIELIIOLOLDLDLDLILOLSOLULSIODLI I eee 
eS 
G wWEITS (4,544) 
C RERITE (5,344) 
G DO 7134 I=1,N 
C WR —— {Ol fe) = ae 


CI (33 Set TE (8, 90)  (e (1 7, ee 


C 

é mex HAS SLED OF FILIZR REACHED EOINI 25? 
Th ( (¥=25) oNEL OC) GOton 222 

ec 
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T4E1=. 50*ATAN( Z*PHKH1(1,2) / (PAK41 (1,1) -PKK 41 (2,2))) 
Gree a ee yee 11 (2, 2) fz. FER < 41 (1,2) 

1 /SIN(2.*TEEF1) 

SiGe Y= (eet (1) +PAx% 112, 2) ) /2.-PKAN1 (1,2) 

ie SUN eee THE 1) 

ee See) ee eS 

Y=q(SIEG2Y) **s 5) <5 

Pi=3.14159265/12 

@1=CCS tiae |) 

Seo rine 1) 

DGmISCMTEELIE=1,25 

el eee 

a lille ke ia dali aia 
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a 


Eas Pon (a, 1682) Xf (2 ELL IP) ,YP (1ELLIP) 
1992 FORMAT (2F15.7) 

C eee END OF ELLIPSE CALCULATICN 

C 


Pee CONTINUE 


C Peete E KM SAUTE 
WPITE(8,22254) 

merc tAr(//," 1X472(J)), TkK443 (Jd), ITKMZ4 (JJ) 
merece. Jy 0) 
Dem) hae) 
yd I= 1,5) 

RITE (8,2225) 

eee FORMAT (/7,* TKP7I2(44), THEGS (4M), TKP24 (4H), 
fe tAcouat) ©) 
APITE(d,2223) (TK272 (MM) , TK P48 (AS), TKP24 (MM) ,TK26 (44) 
P= aN) 
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ek ok ok ee he ok Ato ee ok oo eK ke ok oe ee eK ee ee EO OO Ee 
C Hee (a ee) 
WRITE (EG, 393) (TIMES (1), LAT (1) , LONG (1) eee 
1 Vee 
WRITE(d,393) (BTisE&(1),ELAT (1) , EUONG( Deere 
1, XK 46 (1), L= mz) 
393 FOSMALT 4 1Gpor ace) 
e Were Dd gaa) 
WRITE (4,393) (TIME(1),LAI (1) , LON (1)) vo 
1 pia ty 82) 
WRITE(4,393) (BTIA4E (1) ,ELAT(1) , 21630 ( aes 
Pi Pee el) ga = ie 
c 


ITNTIAS=INT (TIME (4) ) 

RLTIAE=FLCAI (INILAE) 

[fF (RLDIAE. UES Die) ere coe 

CALL TEMR(EBLA1, LONG, TIME, BTIAS, X87 ee, Deed 


339 CONTINUE 


C tee YRITE FOET ALD #082 TC FIER "Ceo 
WEITE(€, 329) 
329 FORAAT(//,5%, JULIAN ER £021 
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fae es we ete et), ec AT (zt), BLCNG (1), =2Y (1) ,2X (4) 
1 p= pee} 
oe es eee ad Ik yn t, 10) 
ie os (et et), 8 = 1, 10) 

9898 FORMAT (F10.4) 

ParTE (oe4 10) 

410 FORMAT (*'1") 


C ALTE (4a, 7777) 
S 
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C eee Weer eGo) ki, 2) 
er seta tal) 
C WRLIE(4,7777) 
c CA timeorere (1 iA, Yh Nein |) 
Cc Cee eeemee at Le Ye iN, S) 
LINN+2 
Bence (hit) = 109 
eco y= 150 
LAgre (Lb 41) = 
MAL (Lis o0 
ARITE(4,7777) 
c Gibitae eee (LONG, LAG, LE, 1) 
C Cibolo Merrie ra, iN, 2) 
C CAhee soe etiOGNG, LAL, iLL ,)) 
c Chie tPA, Y aoe, — ) 
tC 
eer ue ENG C ELCT Filt FOR DISSPLA SXF KEES 
. MUR “op evA LES 
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WRLTE(19,2423) 
£429 FCRMAL (*INMEX',3n,°T142", tee 


S 
*,6X, 'YHR',5%,°LOMe" , 5X, ie) eee 


e 
E VALUES JOM umenne aines 
WRITE (19,2430) (I, TITHE (1),G11(1),G31{1),)¥H 01) ee 
*LCNG (I) ,LAI(I) ,XH(1), I = 1,4) 
C 
2430 ECRMAT (14,2859. 2) 
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ar LE (ere ienienis 
S WRIIE(9, 1016) k 
DO TGs =i. 
ARETE (8,90) (6 (07 oe Fe 
6S CONTINUE 
ARITE (9,10) 
16 FORMAT (' ',7,3%,'CCV. MAI. OF PR 
cc 145 J=1,N 
Prats JO) (ore eee 1 
13 CONTINUE 
Sud sz 
END 
CHAN WAWAN Ad SMM WA WS MW Md WM MMM WB 
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PHI(In,IC) 
TELL (0a, 1c 
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Cope ee) 
DO ie = 
ene, Wee = 
Fo= Fel. 
DO 13 Ik 
DemSMLe 
IF (A2S (CCE 


Bibs iC) 
P2lUaN 


Ee 


Meow Uose le TNE CUuNlePurL Es 


ec vVArITANCS 


SJSsKOULiNe 


T/E*C (15,1C) 
PHI (IR, 1C)+CCR(I8,I<) 
eine, Tey eT /( (ht 1c 


pes ee 2)! ) 
ep 
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1,N 
1,4 
QO. 
es 


eta Cer h (ES ky SCOR (A, 1 C) 


Poeaia bey ti7 (sri. ) *COn(IR,IC) 


(IR,1C)) -GI.TEST*ABS (PHI (Ie,IC))) 


eee, agp kL, i, LE} 
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Do 
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JO 
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TPO (25 meee 
>HI (12,12) ,PHET (12, 12) , eee een 

) = P(K/K-1) *H1¥* (8*P (X/K-1) ET + 2) 
PHI*®P(A-1/K-1) *PHIT + © 

(K. ES. Ty (Ger ners 

CALL TRANS (PHI,N,N,CHII, ND, MD) 
Gal T 


ALL ®ROD (KK, 2 HIT,N,N,N, DESO oP eewney 
CALL PROD (ERI,TPMUP, N,N,N, cece tice een 
CALL ADD (12421,9,N,N, Es sot lee 
ARITE (3,555) 

FORMAT (/! AATRIX EKKM1 °) 
DCMS0ee T= 1 
HarTPh (6,90) (Pawel (aed) = ae 


NEL Sa, Sop es) sien) eae 
GouGe 
CXL Pats th ie pei) 


FORMAT (! il ') 

BC 20 2-<t =a 

WRETE (e790) (Ht) lee 
FORMAT (127E11.4) 

ASITE(8, 36) 

FUSMAT (' HT =f) 


woate (8790) (hi 70) See 

CALL PSGD (2KK41,HT,N,N, Mp LENS, MC Moe LE) 
CALL 220D(H, 1E%P,43,N,M, 1e121,N on 
GRITE(8, 34) 

FORMAT(’ & 2 HI ‘) 
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TEASS(2,2)/DET 
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alee i = 


WRITE (38,31) 
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meee eos, 1) / DEL 
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4) 
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eae (KK) 
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mend Oe ey) ) 


30 FCRMYAT({' 
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2 


2 


ener cl eel 
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WRITE (6, 37) 
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WHERE B(K/k) 
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CALL PROD (1 
RETURN 
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SJuRCUTLN 
DI4ENSION 
DO | 2.8 
DO 2a 
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52 Je (ee0) = 


V1 
Nh 


C(I,J) = 
RETURN 
END 


=e! 1) 


F ADD(A,3,N,M,C,ND,4D) 


4RITE (9,73) (te? oa) 
HZnE PKXM1(I,d) = 
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ieee 
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A(ND, 43), 5 (1S Mb) OO (Noe Ln 


=n 


= 1,H 
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eer?) ae 


B (1,9) 
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SUBROUTINE TRANS (A,N,M,C,N2,MD) 
DIMENSICN A(ND,MD),C(XD,NB) 
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SUPROUTINE CONST (o,A,N,%,C,ND,4D) 
DIMENSICN A(ND,™“C),C(ND,MD) 
eel ty 1) 14 


pO 100 I = 1,N 
po 100 J = 1,4 
Coe 0.0 
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Cc; 
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me Aa 
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yt 
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Cliyd) = =A 
RETIS 

Dost ane 
DO 160 J = 1,# 
C(I,J) = C¥*A(I,d) 
BEIUAN 

END 


SHOR CEI NE RROLSsN, BoA ee 
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Dome me ee= Gel oH 
Meise e oO.) GO TC Jo 
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DO Sey = ee 
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pO 35 J = 1,8 
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CUNTINUE 

CONTINUE 
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Dommiemie= Trp 1, % 
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aml J) = Sf (1 Wyo) 
ZEIURN 

Epa 2 
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RETUBN 

END 


L ie od 
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SUBROUTINE SVs yy ee 
DIMENSICN A CNDENT) Siler es aD 
DO 10: = ave 


WRITE (4,20) {i ,ee l,i, eee 
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END 
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Pec UL (oye DSIUL(I))GOTC 10 
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CONTINUE 
neers 2S 0) 
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